LOWER ORDER TERMS IN THE FULL MOMENT 
CONJECTURE FOR THE RIEMANN ZETA FUNCTION 



J. B. CONREYi'2, D. W. FARMERi, J. P. KEATING^, 
M. O. RUBINSTEIN^ AND N. C. SNAITH^ 

1. AMERICAN INSTITUTE OF MATHEMATICS 
360 PORTAGE AVENUE 
PALO ALTO, GA 94306 
USA 

2. SCHOOL OF MATHEMATICS 
UNIVERSITY OF BRISTOL 
BRISTOL BS8 ITW 
UK 

3. PURE MATHEMATICS 
UNIVERSITY OF WATERLOO 
200 UNIVERSITY AVE W 
WATERLOO, ON, CANADA 
N2L 3G1 



Abstract. We describe an algorithm for obtaining explicit expressions for 
lower terms for the conjectured full asymptotics of the moments of the Riemann 
zeta function, and give two distinct methods for obtaining numerical values of 
these coefficients. We also provide some numerical evidence in favour of the 
conjecture. 



1. Introduction 

In [CFKRS] the authors propose conjectures for the full asymptotics of the mo- 
ments of L-functions. A sample conjecture states, for integer k, that the 2A;th 
moment of |C| on the half line can be estimated using a polynomial Pk of degree 
fc^, with the polynomial given implicity as a 2fc-fold residue (see (1.2) below). 

The leading term in the conjecture agrees with the Keating- Snaith conjecture 
[KS] for the leading asymptotics of the moments of C,. Besides that, all the terms 
of the polynomial obtained agrees with known theorems for A; = 1, 2 [I] [H-B], and 
the leading term agrees with conjectures made through earlier and distinct number 
theoretical methods for fc = 3,4 [CG] [CGo]. 

The method used in [CFKRS] to conjecture the full asymptotics relies on number 
theoretic heuristics based on the approximate functional equation. The conjecture 
is supported by the fact that the formula coincides with an analogous expression 
in random matrix theory [CFKRS2] for the moments of characteristic polynomials 
from the unitary group, the main difference being that the moments of C have extra 
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arithmetic information that does not show up for random matrices. Perhaps the 
most compelhng support for the conjecture, though, are numerics that confirm the 
conjectured moments. 

For those numerical confirmations it is necessary to use all terms arising in our 
heuristics. The purpose of this paper is to investigate the lower order degree terms 
which appear in the conjectured polynomials. Specifically, we 

(1) describe an algorithm to obtain meromorphic expressions in k for the co- 
efficients of the polynomial Pk{x). Our main results are given in Theo- 
rems 1.2-1.4. 

(2) explain how one can numerically compute the lower order terms and to 
provide further experimental confirmation of the full moment conjecture, 
including for non- integer values of k. Numerical values for the coefficients 
of these polynomials for k = 1,2, ... ,7 are listed in [CFKRS] without 
explanation, with some numerics confirming the conjecture for A; = 3. 

At the end of the paper we also outline the analogous approach for moments 
of quadratic Dirichlet L-functions and of quadratic twists of an elliptic curve L- 
function, in both cases evaluated at the critical point. These two cases are examples 
of unitary-symplectic and orthogonal families respectively [KaS] [KeS2] [CF] . 

Before stating our results, we introduce notation and conjectures from [CFKRS]. 



1.1. Moment conjecture for Let 



A{zi,... ,Zm) = Yi 



Zi) 



l<i<j<m 



denote the Vandermonde determinant. 



Conjecture (see [CFKRS]). 

For positive integer k, and any e > 0, 

^ |C(l/2 + it)\'''dt = r Pk (log ^) dt + 0{T^/^+% (1.1) 
Jo 

with the constant in the O term depending on k and e, where P^ is the polynomial 
of degree fc^ given implicitly by the 2fc-fold residue 



Pk{x) = 



(-1)' 



(27ri)2fc 



/ . . . / G{z,,...,z,,)^IS?{z,,...^z,,)^^ - dz, . 

TT .2fc 



n 



with the path of integration over small circles about Zi = 0, where 



k k 

G{zi, ... , Z2k) = Ak{zi, ... , Z2k) n n ''(^ + ^» ~ (1-3) 
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and Ak is the Euler product 

Ak{zi, ... , Z2k) 



P ij=l .7=1 ^ P 



' ' ' e{9) V e{-9) ^ ' 



1 - — 1 de 



(1.4) 



+ k-Zr, 



= TT y TT . (1.5) 

Here e{e) = exp(27ri^). 

We use both these expressions for the local factor for A/.. The first is used in 
obtaining meromorphic expressions in k for the coefficients of Pk{x). 

The second expression, derived in [CFKRS][2.6], is used to numerically compute 
Ak{zi, . . . , Z2k) for specific values of zi, . . . , Z2k- The individual terms in the sum 
over j in (1.5) have poles (though these poles cancel out when summed over j, see 
the paragraph following [CFKRS][2.6.16]) and when we numerically evaluate these 
terms individually, we take care to avoid the poles by making sure that the -Zj+fe's 
are distinct. 

The main point of the conjecture is that we believe it gives the full asymptotics of 
the moments of zeta. While our numerical results in Section 5 are consistent with a 
remainder of size 0(T^/^+'^), there is some debate regarding the error, especially in 
relation to moments of other families of L-functions [CFKRS] [Z] , and it would be 
worthwile to carry out more detailed testing concerning the nature of the remainder. 

The leading coefficient of Pk{x) will be shown in Section 2.1 to equal 

fc-i ., 

ak 

with 



ak 



Hil - p-'f 2F,{k,k;l;l/p) (1.7) 



and 2-Fi(a, b; c;t) the Gauss hypergeometric function 

F(n h-r-f\- ^("^ ^ na + n)T{b + n) 

2l'i{a, b, c, t) - > — — — -. (1.8) 

^W^W^ r(c + n) n! 

This agrees with the leading term that was first conjectured by Keating and Snaith [KeS] . 
More generally, letting 

Zis) = x(s)-^/\(s) 

with 

we conjectured [CFKRS] for shifted moments 

T ^ 
/ Z{l/2 + it + ai)--- Z{l/2 + it + a2k)dt^ [ Pk (a, log ^) dt, 
{ (1-9) 
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where 



i=i j=i 

with the path of integration being small circles surrounding the poles a^, and 
— 1/4 < ^aj. One recovers the moments of ( by setting the shifts equal to 
0, and observing that Z{s) — Z{1 — s). 

An alternative formulation of this conjecture also given in [CFKRS] involves a 
combinatorial sum and is established by the following lemma. 

Lemma 1.1. 5*66 [CFKRS], Section 2.5. 

Suppose F{u; v) = F{ui, . . . , w^; vi, . . . , Vk) is a function of 2k variables, sym- 
metric with respect to the first k variables and also symmetric with respect to the 
second set of k variables. Suppose also that F is regular near (0, . . . , 0), and that 
f{s) has a simple pole of residue 1 at s = but is otherwise analytic in a neigh- 
bourhood about s = 0. Let 

k k 

H{ui, . . . ,Uk;vi, ...vk)= F{ui, ... ,Vk)Y[Y\. ■f^'^^ ~ "^-j)" 

If for all 1 < i,j < k, ai — aj+k contained in the region of analyticity of f{s) 
then 

(-1)" 1 J £ H(z,,.. .z^,)A{z,,....z.,,f , 

E^GH-f^("<T(l):--- :"rx(2A-)), (l-H) 

where one integrates about small circles enclosing the aj 's, and where S is the set of 
{^k) permutations a S S2k such that < • • • < a{k) and a{k+ 1) < • • • < a(2k). 

Equation (1.2) allows us to obtain the coefficients of Pk{x) by computing power 
series expansions and then the residue of the r.h.s, giving meromorphic expressions 
in k for the coefficients which can also be evaluated to high precision numerically, 
even for non integer fc. In practice we have been able to do so for the first ten coef- 
ficients of Pk{x). If A; G Z, fc > 3, to obtain numerical values for all fc^ coefficients of 
Pk{x) we developed a second method using equation (1.11). This involved taking 
small distinct shifts and high working precision to capture cancellation amongst 
the order fc^ poles of the r.h.s. of (1.11). 

1.2. Results. Our first theorem below explictly gives the coefficients of Pk{x) in 
the full moment conjecture for the Riemann zeta function. These are described in 
terms of the multivariate Taylor coefficients of 

—Ak{zi,...,Z2k) TT {zi- Zj+k)C{^ + Zi- Zj+k). (1.12) 

CLk 

l<h3<k 

We let &fe(a; /3) denote the coefficient of z"^ . . . z'^''z^^-^ . . . z^k in the Taylor series 
of (1.12). Let 

fc 



(1-13) 



LOWER ORDER TERMS 



5 



and likewise for /?. Notice that the function in (1.12) is symmetric in zi, . . . , Zk and 
in Zk+i, ... , Z2k, and that 6fe(a; /3) = (— l)l"l+l'^l6fe(/3; a). We may therefore collect 
together terms in the Taylor series accordingly and express (1.12) as 

^ 6fe(a; /?) [zr ■ ■ ■ zl^zlX^ • • • 4 ± sym) (1.14) 

The 'sym' indicates that we group terms that have exponents of the same form. 
Thus the sum over a; /3 follows the convention that a > /3 lexicographically, we list 
the Qi's and /3j's in decreasing order, and we supress the Oj's and /3i's that are 0. 
For example, all the terms of degree 4 are collected with coefficients: 6/c(l, 1, 1, 1; ), 
&fc(2, 1, 1; ), 6fc(2, 2; ), 6^(3, 1; ), 6^(4; ), 6^(1, 1, 1; 1), 6^(2, 1; 1), 6^(3; 1), 6^(1, 1; 1, 1), 
6fc(2; 1, 1), 6fc(2; 2). The terms that go with 6fe(l, 1, 1, 1; ) are 

^ ^ •^12-^^3 -^14 ~t~ -^fe+ii •^fc+i2 -^fe+is •^fe+i4 • 

I<ji<i2<i3<u<fe 

Theorem 1.2. Let Pk{x) be given by equation (1-2). Writing 

Pk{x) = co{k)x''^ +ci{k)x''^~'^ + ... + Ck2{k), (1.15) 

we have 



k— 1 

c.(fc) = afe n T^^nl E 2i-^("'/5)6fe(a;/3)iVfe(a;/3), 

'=0 ^ \a\ + \l3\=r 



(1.16) 



where ak is given by (1-7), the function 6 [a; (3) equals zero unless a = (3 in which 
case it equals 1, and Nk{a;(3) is defined by 

\ /— / rearrangements ( J. X I ) 

a, T of a and /3 \ ' / 

The function Mk is given as a 2k x 2k determinant in equation (2.63). 

By rearrangements, we mean distinct permutations. Two permuatations a {a) 
and /z(a) are said to be distinct if a^; ^ a^. for some i. For example, if ai = 7, a2 = 
5,0:3 = 5 then the two permutations 02,0:1, as and 03,01,0:2 are not distinct and 
would only be counted once in (1.17). 

The reason for writing Nk{a; /3) as we have, with the factor ^nf=o^ {k+i)\ ) ' 
explained by the next theorem. 

Theorem 1.3. Nk{a; f3) is a polynomial in k of degree < 2(\a\ + |/3|). 

This theorem allows us to determine Nk{a;P) explicitly by evaluating (1.17) at 
2(|o| + \/3\) + 1 values of k and interpolating. A few example A^fc(o;/3)'s are given 
in (2.70). 

Finally, the coefficients bk{a; 0) that appear in Theorem 1.2 can also be explicitly 
determined. 

Theorem 1.4. The Taylor coefficients bk{a;P) of (1.12) can be written explic- 
itly as a polynomial in: k, the Taylor coefficients of sC(l + s), and, the Taylor 
coefficients of log{Ak{zi, . . . ,Z2k))- The latter Taylor coefficients can further be 
expressed explicitly as a sum over all primes p of a rational function in: k, p, 
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log(p), and finitely many Gauss hypergeomtric functions 2-Fi(A;i, fc2; to; where 
ki,k2, TO S Z, ki,k2 > k and to > 1. 

To illustrate what the last theorem looks like in practice, see equations (2.51) 
and (2.43). 

This paper is structured as follows. In Section 2.1 we prove Theorems 1.2 and 1.4, 
and also give a procedure to determine the polynomial and rational functions of 
Theorem 1.4. In Section 3 wc prove Theorem 1.3. 

Section 4 is devoted to numerical evaluation of the lower order terms. Two 
different methods are described. The first involves numerically computing the terms 
that appear in in Theorem 1.2, while the second uses (1.11), small shifts, and very 
high precision to capture cancellation amongst the poles of the summand. Data 
supporting the full moment conjecture is then presented in Section 5. 

We also provide plots of the coefficients Cr{k) as a fimction of k for r < 7 and 
also of the zeros of the polynomials Pk{x) for several values of k. 

In Section 6 we briefly describe the analogous approach for quadratic Dirichlet 
i-functions and of quadratic twists of an elliptic curve L-function. 

2. Lower order terms in the moments of ( 

2.1. Evaluating the residue explicitly. The 2fc-fold residue in (1.2) involves ex- 
tracting the coefficient of Hi'" -2^'^"^) i-e. a polynomial of degree 2k{2k — 1). The 

Vandcrmonde determinant has degree 2(^2*^) = 2k{2k — 1). However, the prod- 
uct Yli=i rij^i C(l + Zi — Zj+k) has poles which cancel fc^ of the Vandcrmonde 
factors. Hence, in (1.2), we need only take terms in the Taylor expansion of 
exp ^1 Xli=i ^» ~ ^j+fe) up to degree k'^. Therefore, Pk{x) is indeed a polynomial 
of degree fc^ in x and we write 

Pk{x) = co(fc):c'=' + ci(A;)x'='-i + C2(fc)a;'='-2 + • • • + Cfc2(A;). (2.1) 

One complication in developing expressions in k for the Cr(fc)'s is that the Van- 
dcrmonde determinant in (1.2) prevents us from separating the integrals. However, 
this can be overcome by introducing extra variables and pulling out the Vandcr- 
monde as a differential operator. We illustrate the method for the leading term 
Co(/c) and then generalize. 

Noticing that ^^(0,0,... ,0) = Ofc (set all the variables to in (1.4), apply 
Lemma 2.3, and compare to (1.7)), the leading term is given by 

Co(fc)x''' = ^ {2-KiY^ / ' " / ■ ■ ■ ' • • • ' Zk)i^{Zk+l, ... , Z2k) 

cxp( f ~ Zi+k] 

X ^ 5j Uzi . . . dz2k (2.2) 

i=l 

(the fc^ poles of the ( product have sign opposite from the Vandcrmonde factors 

that they cancel, and these k'^ minuses cancel the ( — 1)'^ in (1.2)). Comparing the 
degree of the Vandcrmonde determinants in the numerator, with the degree of the 
denominator shows that only terms of degree fc^ in the taylor expansion of the exp 
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contribute to the residue. Changing variables Ui — xzi/2 and then relabeling m 
with Zi gives 



X Cfe 1 



j> ■■■ j> A(2;i,... ,Z2fe)A(zi, . . . ,Zfc)A(2;;;+i, . . . ,Z2k) 

n 



— I I. I 

2^ dzi . . . dz2k- (2.3) 

i=i 

Introducing extra variables Xi, we consider 

1 / / exp feiii a^i^i) 
Tifc f ■■■ f pizi,... ,Z2k) ^ '-dzi...dz2k- (2.4) 



{2my 



with p(2;) a polynomial in 2:1, . . . , Z2k- Pulling out the polynomial p, (2.4) equals 



1 f /• cxp 



p{d/dxi, ... , d/dx2k) .„ .^2fc ^:jT ^cJ^i • • • dz2k, 

^^'''> ■> ■> " .fc (2.5) 



2fc 

n 



2. 



and taking the residue gives 



2fe 2fe-l 

Therefore, (2.3) equals 

fe^ 2/c 2k-l 

^^^q{d/dxu... ,d/dx2k)^^^^, (2.7) 
evaluated at a;i = . . . = a;fe = 1, Xk+i = . . . = X2k = —1, with 

q{d/dxi, ... , d/dx2k) = A{d/dxi, ... , d/dx2k)A{d/dxi, ... , d/dxk)A{d/dxk+i, ... , d/dx2k)- 

(2.8) 

Two lemmas allow us to reduce this further. 
Lemma 2.1. 

n 

A{d/dx„... ,d/dXr,)r[f{x,)= f^'-'\xi) 

^ nxn. 

Proof. This follows using the definition of the Vandermonde determinant 
A{d/dxi, d/dxn) = d^-^/dxi'^ 

nxn. 

Noticing that row i of the matrix only involves Xi, we factor the product into the 
determinant. □ 

We can now consider the effect of applying the three Vandermonde's in (2.8). 
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2k 



Lemma 2.2. 

A{d/dxi, . . . ,d/dx2k)A{d/dxi, . . . ,d/dxk)A{d/dxk+i, ■ ■ ■ ,d/dx2k)Y]_f{xi) 

1 

evaluated at xi = . . . = Xk = 1, Xk+i = . . . = X2k = —1 equals 



/(I) 












/(fc-l)(l) 




/(3fe-2)(-2) 


/(-I) 



















2fcx2fe. 

("i/ie /irsf fc rows of this 2k x 2k matrix involve f and its derivatives evaluated at 1, 
while the last k rows have the entries evaluated at —1). 

Proof. Consider first what happens when we apply just the last two A's. By 
Lemma 2.1, 

2k 

A{d/dxi, ... , d/dxk)A{d/dxk+i, ... , d/dx2k) ]J f{xi) 



^kxk Ofcxfc 
Ofexfe -Sfexfe 



(2.9) 



with 



Aij = P'-'\x,) (2.10) 

S.,, = f^'-^\x,+k). (2.11) 

Expanding this determinant, we get a sum of (fc!)^ terms each of which is a product 
of the form 

k 

sgn(/z) sgn(i.) []/(^'-i)(xi)/(''--^H^i+fc) (2.12) 

where Hi, . . . , ^k and vi, . . . ,Vk are permutations of the numbers 1, 2, . . . ,k. Ap- 
plying the third Vandermonde A{d/dxi, . . . , d/dx2k) to a typical such term gives, 
by Lemma 2.1, 



sgn(/i) sgn(i/) 





-'Hxi) 


f^'Kxi) . 


. /(''i+2'=-2)(xi) 




-^\X2) 


&^\X2) . 


. /(''=+2'=-2)(x2) 




-^Kxu) 


&'Kxk) ■ 


f(^,+2k-2)^^^^ 






f^'^^Xk+l) . 


/(^i+2fe-2)^^^^^ 




^\xk+2) 


f^-'Hxk+2) . 


. f(''^+^'^-'\xk+2 




-^\x2k) 


f^-"'\x2k) . 


. /('^^+''=-'n^2fc) 



(2.13) 



2fex2fe. 
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Setting xi = . . . = Xk = Xk+i = ... = X2k = —1, wc may rearrange the 
first k rows and the last k rows so as to undo the permutations fx and v. This 
introduces another sgn(/x) sgn{i/) in front of the the determinant. Hence each such 
term contributes the same amount, 



/(I) 


















/(-I) 






/«(-!) 






/(/=-!) (_1) 




j-(3fc-2)(_;^) 



(2.14) 



2kx2k, 



and summing over the (fc!)^ pairs yu, J/ gives us the Lemma. 
Applying Lemma 2.2 to (2.7) yields 



□ 



co(fc) = ^ 



r(2fc)-i 

r(2A;- l)-i 
r(A: + l)-i 

-r(2fc)-i 
r(2fc- i)-i 



r(2fc- i)-i 
r(2fc-2)-i 

r(fc)-^ 
r(2fc-i)-i 
-r(2fc-2)-i 



(-i)'=^r(fc + i)-i (-i)'=+ir(fc)-i 



r(i)-i 
r(o)-i 

r(-fc + 2)-i 
r(i)-i 
-r(o)-i 

(-i)3'=-ir(-fc + 2)-i 



(2.15) 



2fcx2fc. 



}. The first A; rows and last A; 



Here, we take T{m)-^ = if m e {0, -1, -2, -3, . . . 
rows are identical except for the presence of a checkerboard pattern of minus ones 
in the latter rows. The i,j entry above equals 



r(2A;-i-i + 2)-i, ifl<i<fc; 
(_l)i+j-fc-ir(3A; + 2)-i if fc + 1 < i < 2fc; 



(2.16) 



We show later that (2.15) equals nf=o^ l\/{k + l)\. See Lemmas 3.1 and 3.2, with 

Ci = fi = 0, Ci = k + i - 1. 

Next we consider in (2.1) the rth term of our polynomial P^ix). To evaluate 
Cr{k) we examine the power series expansion of the integrand in (1.2). As in our 
consideration of co(fc), wc first cancel the poles of HiLi 11^=1 CO-+^i~^j+k) against 
the Vandermondc, and write the integral in (1.2) as 



Pkix) = 



1 



1 



(27ri)2*= 



A{zi, ... , Z2ft)A(zi, . . . , Zk)A{zk+i, ... , Z2k) 



2k 

n 



,2fe 



Ak{zi,... ,Z2k) Yi i^i ~ ^3+k)C{^ + - Zj+k) 
l<ij<k 



e2 



^i+kdzi...dZ2k- 



(2.17) 
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Because of the various symmetries satistfied by the factors of the integrand, our 
job of determining the series expansion of 

Ak{zi,... ,Z'2k) n (^i - ^j+fc)C(l + - ^j+fc) (2-18) 

l<i,j<k 

is not as difficult as might be supposed. 

2.1.1. Series for Ak{zi, . . . , Z2k)- First, wo write 

Afe(zi,... ,Z2k) = CXp (log(AA;(zi, . . . ,Z2fe)) (2-19) 

to turn the Euler product defining in (1.5) into a sum over primes. Obtaining the 
series for \og{Ak{zi, . . . , Z2k)) will allow us, in conjunction with the method in Sec- 
tion 2.1.3 for multiplying scries, to recover the scries for cxp (log(Afe(zi, . . . , ^2fe))) 

Because Ak{zi, . . . , Z2k) is symmetric in zi, . . . ,Zk and separately in Zfe+i, . . . , Z2ki 
we distinguish two sets of variables, and let 

Bfe(ai,a2,... ,ak;0i,p2,... ,Pk) (2.20) 
denote the coefBcient of a typical 



Zl 



^2fe 



ai! afc! A! /3fe! ^^'^^^ 

in the multivariate Taylor expansion oi\og{Ak{zi, . . . , Z2k))- Here, we prefer to keep 
factorials in the denominator rather than absorb them into Bk for convenience in 
describing the procedure to obtain the coefficients through differentiation. 

By the above mentioned symmetry, we use the convention, when writing Bk{a; (3) 
of only listing the non-zero a^'s and /3j's, and writing them in decreasing order. Also, 
because 



we have 



where \a\ 



Ak{zi, ■■■ , Z2k) = Ak{-Zk+l, ■■■ , -Z2k, -zi, ■■■ , -Zk+l) 

Bfc(a;/3) = (-l)W+l/3lBfe(/3;a) 
-'^lOii- Therefore 

Ak{zi, ... , Z2k) = 



(2.22) 
(2.23) 



afecxp 



)(sfe(l;) 
+Sfe(l;l) 



k 

E 

1 



E 

l<i,j<k 



Zi - Zi+k + Bk{l, 1;) 

2k 2 

z^Zj+k+Bk{2;)J2^ 



E 

l<i<j<k 



ZiZj -\- Zi-^kZj-\-k 



+-Bfe(l, 1, 1; ) ^ ZiZjZi - Zi+kZj+kZi+k -\ 

l<i<j<l<k 

(2.24) 

The Uk factor comes from the value of the function at the origin Ak{0, . . . ,0) = Uk. 

Next, let I = l{a) denote the number of non-zero a^'s, and m = m{(3) denote 
the number of non-zero /3j's. Since wc arc interested in extracting Cr{k), we only 
need to consider the power series expansion of Ak up to degree r, i.e. I + m < r. 
Since we are assuming in our evaluation of Bk{a; (3) that the aj's and /3i's are in 
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decreasing order, we focus on the first I Zi's (m Zi+k's respectively) and we have, 
together with (1.5), 



/?) = E n ^ n ^ log ^)) 



p i=l 



=1 "^i-y- 



+k 



(2.25) 



2=0 



where 



„i k 



i<i,j<k j=i (2.26) 

Since we are assuming a;+i = . . . = afe = 0, /3to+i = . . . = /3/s = 0, with I + m <r, 
we may as well immediately set Zr+i = . . . = Zk = 0, Zk+r+i = . . . = Z2k = 0. 
Therefore, 

E log(l-ti+^'-^'=+^) = E logll-t^+'^-'^^'+O 

r 

+ - r) (log(l - ti+"') + log(l - ti-^-+'=)) . (2.27) 



By equivalent, wc mean that both expressions have the same series expansion in z 
up to terms involving just zi, . . . , Zr, though not including the constant which, on 
the l.h.s., equals fc^ log(l — t). The main point in doing this reduction is to get rid 
of the k dependence in the summands. 

A symbolic differentiation package (such as Maple) can then be used to compute 



i=i "^^i j=i "^i+k 



(2.28) 



z=0 



applied to the r.h.s. of (2.27) as a rational function in k, log(t), and t. We list the 
terms up to degree 2: 

E log (1 - = log(l - i) - E - '^+>' 



i<i,j<fe 

tlog(t)2 



5Z ^''^i 



+k 



kt\og{tf 



{i-ty 



(2.29) 



The coefficient above of the X]i<j<j<fc ZiZj + Zi+kZj+k term equals zero. 
Next, applying (2.28) to 

log JJ (l - e((9)ti+^^) (l - e(-6l)ti-^'=+^) d6lj 

we end up, by the chain rule, with a rational expression involving partial derivatives 
of the form 

n |i n ^ / n (i - -v^i^^-r (> - < 

i=l <^^i i=l C^i+fc -^0 -1 



(2.30) 



(2.31) 
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Now, 



—1 °° 



m=0 



SO 



The sum above is of the form 



E 



(2.32) 



OC _1 

^(l- ei9)ti+^) " = log(i)^ Yl iem-^rm". (2.33) 

'^^ ^ ^=0 m=0 



(2.34) 



which can be evaluate by applying [wd/dwY to the geometric series 1/(1 — w) = 
Ylm=o '^^^^ '-^'^ ^® expressed either in terms of Stirling numbers of the second 
kind or, alternatively in terms of Eulerian numbers [St]: 



m=0 



Z=0 



1=0 



(2.35) 



(the latter sum is taken to be 1 if c=0). We prefer to use the latter. Thus, (2.33) 
equals 



c-l 



iog(f)^(i - e{e)t-^)-^-' ^(^' + 1)^)*^'+' 



)/2 



(2.36) 



Likewise, 



z=0 



d-l 



(-log(i))''(l - e(-0)t3)-'i-i^i;(d,/)e(-(^ + l)0)t('+i)/2. (2.37) 



Applying this to (2.31) and expanding out, we need to evaluate integrals of the 
form 



(l_ 6(0)^1/2)-'-^^* (l_e(-0)fV2) 



e{ce) de 



(2.38) 



where C € Z, —^di < C <Y^Ci. The integral above can be expressed in terms 
of Gauss' hypergeometric series. 

Lemma 2.3. Let A,B,C ^1, A,B >l, Q <t <l. IfC>0 then 
^ (l - e(0)ti/2) (l - ei-e)t^/^^ e(C0) (i6l 



B + C- 1 
C 



2Fi(AB + C;C+l;i). 



(2.39) 
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IfC<0, then 



(l - e{9)t^/^) ^(l-e(-^)W2) ^ e{Ce)d9 
= tl^l/^ (^^ + l^i' " ^) 2^^i(i3, A + |C| + 1; t). (2.40) 



,ici/2r^+ici-i 

Proo/. Assume C > 0. We can expand (l - e(6')iV2)-^ and (l - e(-6l)iV2)-^ 
using the binomial series: 

(l-eWtV2)-^ = 1 + Ae(^)tV^ + ^^ieie)tyr + ^ ^ ^ (e(^)tV^)3 + . . . 
(l_e(-0)ti/2)-^ ^ 1 + i3e(-g)tV2 + B{B + 1) (^(_^^^i/2)2 ^ B{B + l)(g + 2) (^^_^^^i/2)3 ^ ^ _ 

Multiply these series together. The integral will pull out the coefficient of e{—C6), 
which equals 

c/2 ( B...{B + C-l) B...{B + C) B...{B + C+\) A{A + l) ^^ 
* V C! + (C + 1)! (C + 2)! 2! * 

= i^/^ (^^ +^ " ^) 2i^i(A, B + C; C + 1; t). (2.41) 

The second formula in the lemma can be obtained by conjugating the first and 
interchanging the role of A and B. 

□ 

Using this lemma, we can write out the Taylor series of (2.30) 
log W {l - e{9)t'^+^^y^ (l - e{-e)ti-"'+^y^ de^ = 

log{,F,ik, k; l;t)) + ^F^(^k,k;l;t) ^ " + 

( -{t\og{t))^k\F,{k+l,k + l;2;tf {tlog{t))^'t')^Mk + 2, k + 2;3;t) \ 
\ 2Fi{k,k-l-tY ^ 2Fi{k,k;l;t) J 

l<i<j<k 

• (Hog(t))2fc22Fi(fc + 1, fc + 1; 2; t)2 



2Fi(fc,fc;l;t)2 



tlog{ty2Fi{k + l,k+l;l;t)\ ^ 
2Fi(fc,fc;l;t) 



^ l<i,j<k 

-{tlog{t)fk\F,{k + l,k+l;2;t)^ {t log{t)f 2F,{k + 2, k + 2;3;t) 



2Fi{k,k;l;t)^ 2Fi{k,k-l;t) 

^log(^)%fl(fc + 2,fc + l;2;^) ^^ ^ z[ 

2Fi(fc,fc;l;t) ^^"^^^ 
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Combining the above with (2.29) we have that the first few coefficients Bk in 
(2.24) are given by 

^ k log(p) log{p)k 2f i(fc + 1, fc + 1; 2; 1/p) 
''^ ^ P-'^ P2Fiik,k;l;l/p) 

. ^ ^ / logjpn^F, (fc + 1, fc + 1; 2; l/pf log(p)^ ("+^) i (fc + 2, fc + 2; 3; f /p) ^ 
' ^[ p^2F,ik,k;l;l/pr p^2F,{k,k;l;l/p) 

pi0g(p)2 /log(p)2fc22F^(fc + l,A; + l;2;l/p)2 



p22Fi(fc,fc;l;l/p)2 



iog(j;j-2Fi(/v + l./,:+l:l;l/y;) 



^^(2;) = -E 



P2i^i(A;,/c;l;f/p) 
kplogjpf / log(p)2fc2^fi(fc + l,fc+l;2;l/p)2 
(p-l)2 +1^ p2 2Fi(fc,fc;l;l/p)2 



log(p)' ( 2 ) (fc + 2, fc + 2; 3; 1/p) log(p)2fc i (fc + 2, fc + 1; 2; 1/p) 
p^2Fi{k,k;l;l/p) p2Fi{k,k;l;l/p) 



(2.43) 



2.1.2. Series for the C product. Let 



s C(l + s) = 1 +70S - 7is' + 1^' - + • • • (2-44) 

be the Laurent expansion of sC(l + s) about s = 0, where the 7„'s generalize Euler's 
constant 



7„= V1^_W'"\ (2.45) 

rn^oo ^ fc n + 1 

K=l 



As with the series for Ak, we can here exploit the symmetries satisfied by the 
product 

Yl {zi-Zj+k)C{l + Zi-Zj+k). (2.46) 

l<i,j<k 

We first set Zr+i = . . . = Zk = 0, Zk+r+i = . . . = Z2k = before applying (2.28), so 
that (2.46) is equivalent, in its series expansion up to terms involving just Zi,. . . , Zr, 
to 



n (z, - z,+k)ca Zj+k) HiziCi^ + 2i))'=-'^(-^.+feC(i - ^.+fe))'-^ 

l<i,j<r i=l (2.47) 
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Again, one may use a symbolic differentiation package to evaluate (2.28) applied to 
the above, and thus obtain the coefficients of the multivariate series expansion 



n ~ Zj+k)C{l + Zi- Zj+k) = 1 + Zi - Zi+k 

l<i,j<k 1 

+l'^k'^ ^i^j + Zi+kZj+k + (271 + 7^ - 7^fc^) ^ ZiZj+k 

\l<i<j<k J l<i,j<k 

2k 2 

+ (72^2 - 72/0-271 A:) ^|- + ... (2.48) 



2.1.3. Multiplying series together. Let us be given two multivariate series of the 
form that appears, for example, in (2.24) 

53 Ck{a; /3) (zr • • • ^r^f;i • • • 4 ± sym) • (2-49) 

a;/3 



In (2.24) we pulled out the constant term from the scries, but generally, the above 
can have a non-zero constant term Cfe(;). The 'sym' indicates that we group to- 
gether terms with exponents of the same form as explained in the introduction. 

We can easily obtain the coefficients of their product by examining, for a given 
a; (3 the various pairs 71; Ai, 72; A2 with 71 + 72 = a, and Ai + A2 = 13. Some care 
is needed in carrying this out. While in collecting terms by Ck{oi]P) we use the 
conventions in the above paragraph, 7 and A need not satisfy 7 > A. For example, 
say with A: = 3, a term of the form z'lz2z'^ can arise through multiplication in 

O/l f 71,1 71,2 71,3 72,1 72,2 72, 3\/ -^1,1 -^1,2 Ai 3 A2,l A2,2 A2,3\\ ..i . 

24 ways as [z-^ Z2 z^ z^ z^ Zq )[zi Z2 z^ z^ z^ Zq )) with 71,1 + 
Ai,i = 3,71,2 + Ai,2 = 1,72,2 + A2,2 = 2, and all the others equal to zero. For each 
of these 24 ways, one needs to look up the corresponding coefficients of both series 
by sorting the 7 and A and possibly swapping, using (2.23), so that 7 > A. 

In this manner, we are able, given the series for Ak{z\, . . . ,Z2k) and the series 
in (2.48), to obtain the series for the second line in (2.17) 



Ak{zi, . . . , Z2k) i^i ~ ^3+k)Ci^ + Zi - Zj+k) 

l<i,j<k 

= akJ2 bk{a; P) (z^ ■ ■ ■ ■ • ■ 4 ± sym) (2.50) 

a;/3 



(to determine the multivariate series for Ak from that of log(^fe(zi, . . . , Z2fe)), one 
uses the Taylor series for exp, applied to log(Aft), and the above multiplication 
algorithm). 
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We list the first few terms: 

Ah{zi,... , Z2k) JJ {Zi - Zj+k)C,{'^ + Zi — Zj+k) 
l<i,j<k 

( ^ 

^ 1 
+(7V + Bfe(l,l;) + Bfc(l;)2 + 27fcSfe(l;)) ^ z,z, ^ Zi+kZ^+k 

l<i<j<k 

l<i,j<k 

^ 2k 

+_(^2^2 _ _ 2^^^ ^ 5^(2; ) + Bk{l; f + 2^kBk{l; )) ^ zf 

1 



(2.51) 



2.2. Determining Cr{k). Extracting the terms of degree r from the Taylor expan- 
sion (2.50) and substituting into (2.17) we have 



k^-r 



ak 1 /' /' A(zi,... ,Z2k)^{zi,--- ,^fc)A(zfc+i, . . . ,z2k) 
k\^2mfkf---f 2k 

i=l 

J2 h{a; /?) ■ ■ • 4' 4+1 ■ ■ ■ 4 ± sym) 

|a| + |/3|=r 

ef Ei=i - ^i+kdzi . . . dz2k, (2.52) 



Now, each term in the second line of the integrand with the exponents of the same 
form integrates the same since the integrand is a symmetric function of 2;i , . . . ,Zk 
and of Zk+i, ■ ■ ■ , Z2k, and also because the contribution from z"^ ■ ■ ■ z^''z^^-^ ■ ■ ■ z^^ 
is the same as for z^^ ■ ■ ■ z^'°z'^\i ■ ■ ■ ^"fe ' '^^'^ ^^^'^ changing variables Uj = 
-Zi and using 6fe(a;/3) = (-l)l"l+l'5l6fc(/3; a). 

Let dk{a;j3) denote the number of terms of a given form of exponent. For 
example, (ife(l; ) = 2fc since there are 2k terms in X]i=i ~ Zi+k- Once can write 
down a formula for dk{a;P) as, up to a factor of 2, a multinomial coefficient in 
terms of the multiplicities of the values assumed by the a^'s and /3i's. Let ma{j) 
denote the number of occurences of j amongst ai, . . . ,ak, and likewise for /3. 

Let Ji denote the largest value amongst the Qj's and J2 the largest value amongst 
the /3i's. Let 



(2.53) 
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We have introduced 6{a;P) to take into account that the ±sym in (2.52) collects 
together the terms corresponding to a; (3 and to /?; a. We have 



2l-S{a;P) 



k — ma{0)\ (k — TOa(O) — TOa(l)\ ~ • ■ • — ITla^Jl — 1)^ 

ma{0)J V TOa(l) J \ TOa(2) / V TOa(Jl) y 

k \/k- m0{O)\ fk - m^(0) - m0{l)\ (k- m0{O) ... - m0{J2 - 1)^ 



= 2i-^("=«(A!)2 n n (2-54) 

since there are (^^(q)) ^^Y^ to choose which -Zj's, 1 < i < fc, have exponent 0, then 
i''^ml{i)^) "^^y^ to decide which of the remaining ^j's have exponent 1, etc, and 
likewise for /?. In the simplification to obtain the second line we used J2j=o "^a(i) — 
k, and similarly for /3. 

Therefore, counting the number of terms that are collected for a given bk{a;(3) 
we get 

^ ' |a| + |/3|=r 

A(Z1,... ,2:2fc)A(zi,... ,Zfc)A(zfc+i,... ,Z2/c) ^c«i M^Pl A 
^ ^1 ■ ■ ■ ^fe ^/s+1 ■ ■ ■ ^2fe 

i=l 

xef S»=i ~ ^'+^dzi . ..dz2k, (2.55) 
Pulling out 

A{zi, ... , Z2/c)A(zi, . . . , zk)A{zk+i, ... , Z2k) E j3)z"^ ■ ■ ■ ^^''^f+i - --z 

H + \0\=r 



X 



from the integral, we have, as in our consideration of co{k), 
•^^f^^) = ^fc^::7^^2(9/(9a;i, . . . , d/dx2k)q{d/dxi, ... , d/dx2k) (-gfc- 1)1 

(2.57) 

evaluated at xi = . . . = = 1, Xk+i = . . . = X2k = —1, where q is given by (2.8) 
and 

q2{zr,..- ,z2k)= E &fc(a;/3)4(a;/3)4'^---2r^f;i---4'- 

|a| + |/3|=r (2.58) 

As in the proof of Lemma 2.2, (2.57) equals 

dk 



q2{d/dxi,... ,d/dx2k)'^9{(J','^) (2-59) 



J2k-1 , 



where g{iJL,v) equals (2.13) (with /(x) = x /(2A; — 1)!) and is comprised of a 
sign and a determinant. The sum is over all k^? pairs of permutations yu, v. 



2k 



(2.56) 
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Applying q2{d/dxi, . . . , d/dx2k) to these determinants we get 



2k^ 



E 

\a\ + \l3\=r IJ., 



^ bk{a, (3)dk{a] 13) / n / m 
^12 sgn(;u) sgn(i/)Mfe(/x, y, a, /3) 



with 



Mk{fJ.,i',a,P) = 

/(^i-i+°i)(a;i) 



(2.60) 



^(Ml+2fe-2+ai)(2;i) 
y(M2+2fc-2+a2)(2.2) 

j:{fj.k+2k-2+ak) (^^^-j 



i'fc+2fc-2+/3fc 



(2.61) 



2fex2fe. 



Setting a;i = . . . = Xfc = 1, Xk+i = . . . = 0:2*; = —1, and rearranging rows (to undo 
the n and v) we get 



2k- 



|a| + |/3|=r o-,r 



with 



Mk{cT{a),T{(3)) = 

r(2fc-l-a<,,)-i 



(-i)E/^- 



r(fc + l-a,J-i 

-r(2fc-/3,j-i 
r(2fc-i-/3,j-i 



r(2A: - 1 - a<,J-i 
r(2fc-2-a<,J-i 

r(fc - a^J"^ 

r(2fc - 1 - /3,j-i 
-r(2fc- 2 -/?,,)-! 



(_i)fcr(fc + 1 - (-i)'=+ir(fc - 



r(i -a,J-i 
r(-aa.)-^ 

r(2-A:-a,J-i 

r(i-/3,j-i 
-r(-/3rj-^ 



(_l)3fe-ir(2-A;-/J^ 



The extra factor of (— 1)^ comes from the extra powers of —1 that arc pulled out 
of the bottom k rows of the matrix. Notice that in order for Mk to be non-zero, the 
ctcTi + i — Vs must be a distinct subset of {0, 1, 2, 2k — 1}, and similarly for the 
P^.+i — I's. The implication of this latter point is discussed further in Lemma 2.4 
below. 

For any a,/3, many of the fc!^ pairs of permutations u, r will give the same 
determinant because of multiplicity amongst the a^'s and /3i's. In fact, since r is 
fixed, most of the a^'s and Pi's will equal zero. As before, let ma{j) denote the 
number of occurrences of j in ai, . . . , afe, and similarly for (3. 

Then (2.62) becomes 

\cc\ + \B\=r ' 3=0 j=0 rearrangements (2.64) 

' ' ' o", T of a and (3 ^ ' 
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By rearrangements, we mean distinct permutations as explained following (1.17). 
Recall that Ji and J2 denote the largest value amongst the a^'s and /Jj's respectively. 
Notice that this simplifies since, by (2.54), 

n -"(j-)! fl -/^O')' = 2'"'^"'^ (2-65) 

is constant, where S{a;(3) is given by (2.53). 
Therefore 

^r{k) = ^ E 2i-*(«=^)6,(a;/3) ^ M,(a(a), r(/3)). 

|a|-|-|/3|=j- rearrangements (2.66) 

' ' ' ' CT, T of a and P ^ ^ 

Expression (2.64) can be pared down further by realizing that 'all of the action' 
takes place in rows k — \a\ + 1, . . . , k and 2k — \P\ + 1, . . . , 2k. By this we mean 
that we need only focus on the rearrangements that have a^i = . . . = Qct^.i^i = 0, 
and = . . . = (i<Tk-\i)\ = 0, because otherwise the determinant will equal zero. 

Lemma 2.4. The determinant Mk{(j{a),T{P)) in (2.63) equals zero, unless a^^ = 
. . . = Oio-j, = /3t-i = . . . = Prk-ipi = which case it might, or might not, equal 
zero). 

Proof. Assume that Mfe((T(a), t(/3)) ^ 0. Let := > 1 be equal to the first 
non zero a^. But this forces 62 ■= otat+s^ to also be > 1, otherwise rows i and 
i + 61 would coincide and the determinant would be zero. But then 83 := a(y^_^_g^_^_^^ 
must also be > 1 otherwise rows i + 5i and i + S\ + 62 would coincide. Continue 
in this fashion until reaching beyond the A:th row, i + Si + S2 + ■ ■ ■ + Sj > k. But, 
Si + ■ ■ ■ Sj < X]m=i ~ \'^\^ SO i + \a\ > k i.e. i > k — \a\. Thus, cto-i = . . . = 
C((7k-\c\ — 0- Similarly, = ... = /3tj._|^| = 0, the only difference in the proof 
being that the rows would coincide up to a factor of ±1. 

□ 

The above lemma greatly improves the speed with which we can evaluate (2.64) 
since all but 0^(1) of the terms can be discarded. Consider 

^ J2 Mfe(a(a),T(/3)). (2.67) 

rearrangements 
cr, X of a and /3 

We will prove in Section 3 using the theory of factorial Schur functions that the 
above is equal to nf=o^ l\/{k + l)\ times a polynomial in k of degree < 2{\a\ + |/3|), 
or else is the polynomial. 

Hence, if we let Ni.{a; /3) denote the polynomial 

CT, T of CK and (3 ^ ' 

we have thus arrived at: 

k—i 

(^r{k)=akYljr^ E 2i-^("'/5)6fe(a;/3)iVfe(a;/3). (2.69) 

The factor Y^=o 1^ {k + l)\is equal to the leading coefficient co(A;). The hk{a, /3)'s 
are the Taylor coefficients of ^Ak{zi, ... , Z2k) Y{i<i,j<k{^i - Zj+k)C{'^ + Zi- zj+k), 
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and the first few arc listed in (2.51). The function S{a; (3) equals zero unless a = P 
in which case it equals 1. We have thus managed to express Cr(fc) as equal to co(fc) 
times a polynomial in k with coefficients linear in the 6fe(a,/3)'s. 

Knowing that Nk{a; f3) is a polynomial of degree < 2(|a| + \f3\) allows us to 
determine it for a given a; f3 by evaluating (2.68) at 2(|q!| + 1/3|) + 1 different values of 
k and writing the unique polynomial of degree < 2(|a| + |/3|) that interpolates those 
values. Since the arithmetic just involves rational numbers it can be performed 
exactly. When evaluating the r.h.s. of (2.68) one should make sure to exploit 
Lemma 2.4 so as to only evaluate Or(l) of the rearrangements. 

In this way, one can find, for example. 



This allows us to write down formulae for ci(fc) and C2{k): 
ci(fc) = Un(^]2fc%(l;) 



and, after simplifying, 



iVfc(l;) 
iVfc(2;) 
A^fe(i,i;) 
A^fc(i;i) 





A;2(fc-l)(fc + l)/2 
-fc^(A;-l)(fc + l). 



(2.70) 



C2ik) = \ak^^J^\k^k-l){k + l){bk{l,l■,)-bk{l■,l)) 
X {2{Bk{l; ) + Jkf - 7^ - 271 + 5^(1, 1; ) - Bfe(l; 1)). 

(2.72) 

The Sfe's above arc given in (2.43). In practice, we were thus able to explicitly 
determine the first nine lower order terms, r < 9. 

3. Proof that Nk{a;l3) is A polynomial. 

Throughout this section we use the following notation. Let ei, 62+1, . . . , ek+k—1 
be distinct integers and /i, /2 + 1, . . . , fk + k — 1 be distinct integers. If /i, /2 + 
1, . . . , /fe + fc — 1 is a subset of 0, 1, . . . , 2fc — 1, let ci, . . . , C/t be the complementary 
subset. Later we will introduce some extra assumptions on the Cj's and /i's, namely 
that most of them are equal to zero. 

The following lemma expresses the kind of 2k x 2k determinant that appears in 
the formula for k determinant involving binomial coefficients. 
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Lemma 3.1. Let Ci and fi be given as above. Then 

r(2fc-ei)-i r(2A; - 1 - ei)-i 

r(2fc-l-e2)-^ r(2fc-2-e2)-^ 



T{k + l-ek)-' 
-V{2k-h)-^ 

r{2k - 1 - /2)-i 



r(fc-efc)-i 
r(2fc - 1 - 
-r(2fc - 2 - /2)-i 



r(i-ei)-i 
r(-e2)-i 



r(2-fc-efe)- 

r(i-/i)-i 



(-i)'=r(fc + i-/fc)-i (-i)fc+ir(fc-/fc 



-r(-/2 



(_i)3/c-ir(2_fc_/,)-i 



+0!(/, + i+0! 

n(fe+0! 



2C!+i-e,+i- 



J Vei+i-1/ 



fcxfe 



2fex2fc. 

(3.1) 



where sgn{f) is (—1) raised to the number of transpositions needed to get /i + 1, /2 + 
2, . . . , fk + k sorted into increasing numerical order. 

Also note, if Ci 's and fi 's are not distinct as required, then two of the rows on 
the l.h.s. will coincide up to sign, and hence the determinant will equal zero. 

Proof. Introducing a {2k — j)l in column j, 1/(6^ + 2 — 1)! in row i and {—iy~^/{fi + 
i — 1)! in row k + i, i = I, . . . , k, the l.h.s. of (3.1) equals 



(-1) 



(''-^1'' \4 ((7 + 1 + /)!(//H 



n 



0! 



l\{k + l)\ 



\e2 + l) 



\ek+k-l) 

_(2k-l\ 
\f2+l) 



(2k-2\ 
\e2+l) 



/ 2k-2 
Vcfc+fe- 



l) 



2fc-2\ 
(2k-2\ 

l/2 + li 



/ 2fc-l \ ,( 2k-2 \ 
\jk+k-l) +\f^+k-l) 

Multiplying the matrix in (3.2) on the right by a unit, i.e. by 

= (-1)^ 



' 0^ 
\ei/ 



{ ° ) 
Ve2 + 1/ 

\ek+k-l) 

+ (A+fc-i) 



(3.2) 



2fex2fe. 



2^ ~ V 2fcx2fe 

the ij entry is, for 1 < i < fc, 1 < j < 2fc 

2fe-l 

E 



m=0 



m 

Ci + i- l 



J-1 

m 



2^- 



i-1 

Bi + i-l 



and the {i + k)j entry is, for 1 < i < fc, 1 < j < 2fc 

2fc-l . X ,x f(-l)/* + 



E 

m=0 



/i + i - 1 



J-1 

m 



i-l 



(-1)^ 







iij^fi + i 
otherwise. 



(3.3) 



(3.4) 



(3.5) 



These binomial identities can be proven by noticing that ('^) (^) = (^) (^Z^) and 
using the binomial theorem. 

Expanding the determinant of this new matrix along the last k rows, and pulling 
out powers of 2's from the first k rows gives the lemma. The sgn(/) that appears in 
the lemma can be obtained as follows. The last k rows of the matrix given by (3.4) 
and (3.5) consist of ±l's with a (— 1)/*+'"^ appearing in row k+i and column fi+i, 
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1 < i < k. We can swap the columns and rows of this matrix so that the bottom 
left k X k submatrix becomes diagonal, and the lower right submatrix becomes, 
Ofcxfe, the zero submatrix. 

There are many ways to do so, but to end up with the determinant in the lemma, 
one should make sure that we do not rearrange the relative ordering of the columns 
corresponding to Ci, . . . , Cfc. If the quantities fi+i appear in increasing numerical 
value, one can simply swap column fi + i with its neighbouring columns on the left, 

one at a time, until it sits in the ith. column. This introduces a (— l)Si=i-'^* into 
the determinant. 

However, if the fi + i's appear out of order, in order to preserve the ordering of 
the columns corresponding to Ci, . . . ,Cfc, one should first swap rows so as to put 
the lower k x 2k submatrix into reduced row echelon form. For example if one has 
fii +ii > fi2 +*2, but ii < 12 then one should swap rows k + ii and k + i2- This has 
the effect of placing the (-l)/n+»i-i -^^ entry {12, fi-, = (12, (fh +n -«2) +^2), 
and the (— l)-''»2+'2-i [j^Qj^ij-y (ii,/i2+*2) = (/j2 +^2 ^ *i) The horizontal 
displacement then needed to get these entries into the Z2nd and iist columns is 
therefore unchanged and equal to fi^ + ft^ ■ Relabeling and repeating if necessary, 
one sees that the contribution to the determinant from the row and and column 
swaps that get the lower left k x k submatrix into diagonal form is 

(_l)sgn(/)+X;/.^ (3.6) 

where sgn(/) accounts for the number of transpositions needed to get /i+l, . . . , fk+ 
k into increasing numerical order. 

One now easily evaluates the determinant by expanding along the lower diagonal 
matrix. This submatrix begins at entry A:+ 1. 1 and this contributes a (—1)'^ to the 
determinant. One also needs to multiply the diagonal entries themselves, and this 
contributes a (-1)^-''*+'"^. 

Collecting the powers of —1 that appear in (3.2) and (3.3) and multiplying by 
(3.6) and by the two factors in the previous paragraph we obtain the sign that 
appears in the lemma. □ 



Now (3.1) equals 

(-l)^(^-)/^sgn(/)n(|^2-.^--^-|(cO....-,L,, (3.7) 

with (z)^ = z{z — 1) . . . (z — /i + 1) the descending factorial, and Xj = ek-j+i- The 
determinant above is essentially a factorial Schur function. 

We have introduced the Aj's and taken the tranpose so as to conveniently apply 
theorems of MacDonald and Chcn-Louck (see [CL] [Theorems 3.2,3.3]) concerning 
the factorial Schur function. The extra (— l)'=('=~i)/2 above comes from swapping 
the i, j entry with the i,k — j + 1 entry. 

In our application to Nk{a-, /?), we found, in Lemma 2.4, that most terms in (2.68) 
can be discarded, and only terms with a^^ = . . . = cc,Tk-\a\ ~ I^ti = • • • = Pt^-ipi = 
contribute. 

So assume that, for some s > 0, ei = . . . = Ck-s = 0, and hence Xg+i = . . . = 
Afe = 0. Next, write Cj = k + j — 1 — cj, where Cj is a non-negative integer. Assume 
that, for some t > 0, ej = for all j > t, i.e. that cj = k + j — 1, if j > t. 
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To apply their theorems, one must first assume that the Aj's are decreasing 
Ai > A2 . . . . We can assume this condition by rearranging the first s columns of 
the matrix in (3.7) if necessary. This will change the sign of the determinant by 
a power of (—1) that depends only on Ai, . . . , Ag. However, since we are actually 
permuting the Xj + k — j's rather than the A^'s, some care is needed. 

Say Xj < Xj+i . We are assuming that the Xj + /c — j 's are distinct (we've assumed 
the Cj + j — I's to be distinct), thus Xj < Xj+i actually implies that 

A, < A,+i - 1 (3.8) 

since otherwise one would have two neighboring Xj + k — j's that were equal. 

Swapping columns j and j + 1, the subscript for the jth column is then Aj+i + 
k — j — 1 = (Aj+i — 1) + k — j, and for the (j + l)st column is then Xj + k — 
j = {Xj + 1) + k — j — 1. Therefore we have replaced (Ai, . . . , A^-, A^+i, . . . ) with 
(Ai; . . . , Aj+i — 1, Aj + 1, . . . ) in which Aj+i — 1 > Aj + 1. Also notice that swapping 
the two columns only permutes the subscripts which therefore remain distinct. 

Continuing in this fashion, we end up with Ai > . . . > A^, with the Aj obtained 
from the Aj's by the above swapping procedure. Notice that 

E^'=E^' (3.9) 
since each swap adds one and subtracts one from the Aj's. 

Theorem 3.2. Assume that Cj = k + j — 1 — Cj, with Cj = if j > t, that 
Ai > A2 . . . > As and that Ag+i = . . . = Afe = 0. Then 

{ci)x.+k-j = {-'^)^^^~^''^'^^{c)x (polynomial in k of degree < 2 Aj j 
Proof. By [CL] [Theorems 3. 2, 3. 3, and page 4150] one has 



kxk 



with c + j - 1 = (ci + j - 1, . . . , Cfc + j - 1), A(c) = lli<i<j<,ni^j - Ci). and 

Wm{z) = E {yi2 -!)••• iVim - TO + 1) (3.12) 

with yi = Zi—i+1 and the conventions that 'Wo{z) = 1, and Wm{z) = if to < 0. The 
main point is that our kxk determinant has been replaced by an s x s determinant. 

Furthermore, 

'Wm{c + j-l)= E (^+i-l-en)(^+j-2-ei2)...(fc+j-m-ei„). 

ii<...<im<fe (3.13) 

The terms in this sum with all i's > t contribute 

^ ]{k+j-l)...{k + j-m) 
m J 

since Ci = ii i > t. This is a polynomial of degree 2m in k. 
The terms with ii < t and the other i's > t contribute 

(m - 1) + - 2) • • • (fc + J - to) + J - 1 - en)' 

^ ' ii=l 
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and the terms with ii < 22 < t contribute 



J^{k+j-3)...{k+j-m) (fc + J-l-e«)(fc + i-2-e,J, 

^ l<ii<i2<t 

both of which are polynomials in k of degree 2m — 1 and 2m — 2 respectively. In 
this fashion one sees that the entries of 

are polynomials in k. Furthermore, expanding this determinant we get a sum of 
products of entries, one from each row and column. However, Wmiz) = if m < 0, 
so only the terms with all Aj — i + j > contribute, and the degree of such a term 
is then 2 ^ A^. □ 

Applying Theorem 3.2 we find that (3.7) equals 



fe-i 



^(c) TT ''f /t^ + 2^'+!-^'+!-' X (polynomial in k of degree < 2y A^). 

to ^•(^ + ')' (3.14) 

Here we have absorbed the factor of ±1 into the polynomial and have also used 
(3.9). 

Thus, given that a^^-^ , Ofo-a + lj • • • , C(ak~^k—1 are distinct and Pn , /3r2+l, • ■ • ) /3rfc + 
fc — 1 are distinct, we get that a typical Mfc(cr(a),T(/3)) appearing in (2.68) equals: 

M,(a(a),r(/3)) = A{^) J] 

X (polynomial in k of degree < 2\a\) (3.15) 
where 7 is the complementary subset of the Pn + i — I's. 

{71,- •• ,lk} = {0,-i-,... ,2k-l}-{Prr,Pr, + l,... ,/3,,+fc-l} 



(3.16) 



and A(7) = UtKjilj - 7i)- 

Now, by Lemma 2.4, most of the tto-^'s and Pn's are 0, 



Q!ai = . . . = a<7fc_|„| = = • • • = Pr^-ini = 0, (3.17) 

and 7j = k + j — 1 — ej where ej = if j > The latter can be seen by noticing 
that + i - I < \P\ + k - 1 so that {fc + |/3|, ... , 2fc - 1} is a subset of 7. Thus, 
starting from the end, Cfe = 0, hence Ck-i =0, ... e\/3\+i = 0. 
Hence, 

A(7)= n ij-i+ei-ej)= n 0'-^) n (i+^^y)- 

l<i<j<fe l<i<j<k l<i<j<k (3.18) 

The first product gives 

fc 

(3.19) 
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and, because — €j — ii i, j > \f3\ the second factor equals 

n (1+^) n (3-20) 

l<i<j< /3 i<i<l/3| 

~ l/3|<J<fc 

However, the product over 1 < z < j < |/3| is a rational number. Furthermore, 
most of the numerator of the product over I < i < \P\,\f3\ < j < k cancels with the 
denominator leaving a polynomial in k of degree |e|. 

But |/3| = |e| which is easily verified as follows. The union of the /3j + i — I's and 
7i's give {0, 1, . . . , 2fc - 1}, so 

k 2k-l 

^A+i-l + 7i= ^ (3.21) 

Substituting ji~k + i — 1— et and simplifying gives A = e,. 
Collecting the above together gives 
k 

A(7) = n ^' ^ (polynomial in k of degree |/3|) (3.22) 

1=1 

(we take the polyomial to be 1 if |/3| = 0). 

Next, we determine the power of 2 appearing in (3.15): 

fc fc 
^li+i - a^i+i -l = ^k + l- ei+i - a^i^^ - l = k'^ + \f3\- \a\. 
1=0 1=0 (3.23) 

Finally, 

{Pti+i + )■ _ pQjyjjQjj^jg^j ^ q£ (jegpgg 1^1^ (3.24) 

1=0 

because = . . . = (^ri^^if^i = 0, and whore wc regard the l.h.s as a fimction of k 
with Prk-j fixed for < j < Therefore, most of the numerator cancels with the 
denominator except for |/3| factors each of which is a polynomial of degree 1 in k. 
We have therefore shown that 

Mk{a{a),T{l3)) = 2'=' TT x (polynomial in fc of degree < 2(|a| + |/3|)). 

f=o(^ + ^)' (3.25) 

Here we have absorbed the extra 2l''l~l"' from (3.23) into the polynomial. This 
proves that Ni.{a;(3) given by (2.68) is a polynomial in k of degree < 2(|a| + \(3\). 

4. Numerical evaluation of Cr(fc) 

Two methods were developed to numerically compute the coefficients Cr{k) of 
the lower order terms. The first relied on (2.69) and we used Maple [M] to take 
advantage of its symbolic capabilities. This approach had the advantage of allowing 
us to obtain the coefficients to many digits precision, and also to make sense of the 
conjecture for non-integer values of fc. This method suffered the disadvantage of 
being difficult to implement, even using a high level symbolic package, and required 
much computational power, so that we only determined Cr{k) in this way up to 
r < 9. This sufficed to compute all the lower terms, for fc = 3, since Pk{x) = 
co{k)x'^ + ci (fc)a;'^ ~^ + . . . + Cfe2 (fc) is a polynomial in x of degree fc^. 
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The second method was comparatively easy to implement, and allowed us to 
obtain many more coefficients. However, it is limited to integer values of k, and 
also presents more difficulties in acceleration therefore yielding lower precision. 

4.1. Method 1. A table of the polynomials Nk{a; P) of degree < 2(|a| + was 
prepared by evaluating (2.68) at slightly more than 2(|a| + |/3|) + 1 values of k and 
interpolating the unique polynomial of said degree fitting those points. The extra 
points were thrown in for good measure as a check against errors. This was done 
for all < |a| + < 9. 

Next a corresponding table of the coefficients bk{a; (3) was prepared, expressed 
symbolically as a polynomial in the jj's and Bk{a; /3)'s as described in Section 2.1, 
with Bk{a', 0) given as a sum over primes, with the summand equal to log(p)l"l"'"l''l 
times a function rational in p and in Gauss hypergeometric functions of the form 
2Fi{k + A,k + B;C; l/p), where A, B, C are non-negative integers, C > 1. A few 
example i?j.(a;/3)'s are listed in (2.43). 

We were then able to obtain, for a given k, numerical values of the coefficients 
Cr{k), for < r < 9. Because Bk{a; (3) is expressed as an infinite sum over primes, 
we used standard methods to accelerate its convergence. Namely, we evaluated the 
first few terms, p < P, to high precision. Then, to evaluate the tail end of the sum, 
p > P, we used Maple's series routine to determined the first few terms of the series 
expansion in 1 /p of the summand, writing it in the form 

where r = |a| + and the dj's depend on the summand hence on a and /3. 
To evaluate a sum of the form 

log(p)'' 



E 



pj 

p>p ^ 

we first wrote it as a full sum minus the front end: 

log(p)'' v;^ log(p) 



(4.2) 



p ^ p<p ^ 

The second sum was evaluated by summing the terms p < P, while first sum was 
computed using Mobius inversion: 



so 



and hence 



00 

i°gC(a) = E-Ei;;^. ^(«)>i. (4-4) 

m=l p 

E4 = E^i°gc(-.), (4.5) 



p" m 

p rn=l 



E^ = (-irf:^(l»SC(m.))«. (4.6) 

p m—1 

Now (log C(tos))'^^^ decreases exponentially fast in m, as can be seen by considering 
its Dirichlet series which is dominated by the first term, and only a handful of m 
on the r.h.s. of (4.6) are needed to evaluate '°^^f^ to a given precision. 
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The factor ak given by (1.7) can be evaluated to high precision in a similar way. 

In this manner we were able to compute Cr{k), < r < 9 for various k. For 
example, for fc = 3 we obtained the coefRcients to about 30 decimal places. The 
actual precision can be predicited from the size of P, as the overall error in us- 
ing (4.1) to approximate the summand for the terms p > P is 0(log(P)''~^/P^). 
In practice, we took larger and larger values of P until the numerics stabilized to a 
precision that we found satisfying. 



4.2. Method 2. The second method we developed to compute c,.(fc) used the com- 
binatorial sum (1.11), small shifts, and very high precision to capture cancellation 
amongst the high order poles of the terms in the sum. Because this method requires 
very little symbolically, this was implemented in C++ using NTL [S] to carry out 
multiprecision arithmetic. 

The basic idea is as follows. The polynomial Pk{x) given by (1.2) can be regarded 
as a special case of the function Pk{a,x) given by (1.10), namely with ai = . . . = 
Q!2fe = 0. One can then use (1.11) to evaluate Pfc(a, x). However, the terms in (1.11) 
have poles if the ai's are not distinct. So, we cannot simply substitute a = 
and sum the terms numerically. Instead we take the limit as a ^ with the 
condition that the ai's are distinct. One must also use very high precision to 
capture cancellation amongst the terms which individually become very large when 
a is small. 

More precisely, let 

(fc \ k k 

Zj+k Ak{zi, ... , Z2k) n n + ~ 
1 / i=ij=i (4.7) 

and let ej = je, where e G C. Then, by (1.11) 

Pk{x) = lim^7?(e^(i),... ,e„^2ky,x), (4.8) 



where ^ is the set of (^^^ ) permutations a G S2k such that cr{l) < ■ ■ ■ < (7{k) and 
a{k + 1) < ■ ■ ■ < a{2k). 

Therefore, expanding exp in its Taylor series, and pulling out the coefficient of 
x'^ we get 

^ 2fe^-'-(fc^-r)! e'^o^^'-^'"W'-- - ^"^-^^ 



where 



k k 



Ck \ ^ ~^ 

y^.Zj - Zj+k ) Ak{zi, ... , Z2k) H H C(l + -Zj - Zj+k)- 

1 / i=li=l (4.10) 

The only complication in evaluating the above for a given k and e is that Ak{zi, . . . , Z2k) 
is expressed as an infinite product over primes (1.4). To evaluate that product, we 
broke it up into p < P and p > P, with P large. For the primes p < P, we 
used (1.5) to evaluate the contribution from p, each factor only requiring finitely 
many arithmetic steps. 
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For the contribution from the larger primes, p > P,wc used a quadratic approx- 
imation for the local factor appearing in (1.4): 

1_ ^ p-2-Zii-2«2+2fe+31+2fc+J2. (4-11) 

l<«l<i2<fc 
1<31<32<*: 

This approximation can be obtained by substituting Uj = p~^/2-^3 and Wj = 
p~^^'^+^''+j in the local factor of (1.4), and working out the terms up to degree four. 
Only terms of even degree appear because the intergral over pulls out just the 
terms with the same number of u's and w's. So, we expand each geometric series 
appearing in the integral over 9 up to degree two, multiply them out, and collect 
terms with the same number of u's and w's. 

The first factor J|(l — UiWj) appears precisely to cancel the terms of degree two 
in the second factor, so wc need only determine the terms of degree four. Noticing 
that the local factor is symmetric separately in the u's and w's, and also if we 
swap u and w, we can determine the terms of degree four by simply computing all 
representative fourth order partial derivatives that have the same number of u's as 
w's, evaluated at Uj = Wj = 0, 1 < j < k. For instance, it is enough to immediately 
set Uj = Wj = if 3 < j < k and then take the partial derivatives: 3^7 gfj af;^ > 
S«rou7^^' 'E^'S^'' evaluated at u\ = ui = w\ = wi = 0. Doing so gives that the 
local factor, up to terms of degree four, equals 

1- X] Ui^Ui^Wj^Wj^, (4.12) 



l<il<i2<fc 
1<J1<J2<'= 



thus giving (4.11). 
Therefore, we used 



n n (1 

1<31<32<'= 

n l<il<i2<'« C(2 + + Zi^ — Zk+jx ^ 2fe+i2)~"^ 
l<jl<J2<'' 

n„<pni<n<'2<^ (1 -p-2-^n-^.2+^'.+«+^'=+.2) 



(4.13) 



to approximate the contribution to (1.4) for the primes p > P. 

To compute (4.9) to D decimal places, we should take e roughly of size 10~^ 
and then use about (r + 1)D digits working precision to account for cancellation 
in (4.9) amongst the order r pole in e of the summands. Since r can be as large as 
fc^, we used [k"^ + 8)D digits working precision, the +8 taken for extra leeway, and 
also chose D to be slightly larger than the desired final precision. 



5. Verifying the full moment conjecture 

In [CFKRS] we presented numerical data supporting the conjecture described 
in Section 1.1 for k = 3. Here we give some more data supporting the conjecture, 
for integer k < 7, and also for several real and complex values of k. Our data 
supports the conjecture, but is not too extensive as our main effort was put towards 
developing ways to evaluate the lower terms rather than to large scale verification of 
the conjecture. Nonetheless, even moderate data strongly supports the conjecture. 
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One experiment we carried out involved comparing the two quantitites 



m + it)\"'dt 



(5.1) 



c 



and 



[ Pk{\og{t/2%))dt, 



(5.2) 



Jc 



for seventeen intervals [C,D\ of length 50000, and k = 3,4,5,6, and 7. We also 
examined the conjecture for several non-integer values of fc, in the latter case inter- 
preting Pk{x) as an asymptotic scries rather than as polynomial. 

For both integer and non-integer k, we used method 1 of Section 4.1 to compute 
the first few lower terms to high precision. For fc = 3, 4, 5, 6, 7, 8 we also computed 
Cr{k) for r < k^ using method 2. 

Tables of the coefficients Cr{k), A: = 3, 4, 5, 6, 7 can be found in [CFKRS]. As k 
increases, it seems from numerics that the first few leading order terms have much 
smaller coefficients than the later lower order terms. For example, when = 4, the 
leading term as listed in [CFKRS] is Co(4) = .24650183919342276 x IQ-^^, compared 
to the largest value Ci3(4) = 38.203306. To verify the full moment conjecture, we 
needed to evaluate Pk{x) at x = log(t/(27r)). Therefore, even if t is moderately 
large, say 10^, the main contribution actually comes from substantially lower order 
terms. In the range we examined, the main contribution for fc = 3, 4, 5, 6, 7 came, 
respectively, from the the terms r = 3, 6, 12, 20, 30 and their immediate neighbours. 

To compare to actual moment data for C, we used Mathematica [Ma] to numeri- 
cally integrate powers of |C(l/2 + it)\ for each interval. Due the oscillatory nature 
of C, we performed the integration between consecutive zeros of C, on the halfline 
using a table of zeros computed with the L-function calculator [R] . 

Tables 1 and 2 gives the values of (5.2) and (5.1) for [C, D] = [50000n, 50000(n-|- 
1)], n = 0, 1, . . . , 16, and k = 3,4, 5, 6, 7. The data for = 3 is a subset of the data 
given in [CFKRS], but otherwise, the data here is new. We see that the pairs of 
columns track one another nicely. 

Figure 1 depicts the difference between (5.2) and (5.1) divided by (5.1) for the 
values in Tables 1 and 2. 

We also present some data for non-integer k, specifically for k = .5, 1.8,3.2, and 
.5 + i taking just the first first few terms, r < 7, of Pk{x). For non- integer k 
we believe, based on our numerics, that Pk{x), no longer a polynomial in x but 
an infinite series, gives an asymptotic expansion for the 2A;th moment of ^, so 
that taking more terms does not necessarily give an improvement. We therefore 
compared (5.1) to 



for a few values of D. We present our data in Tables 3-6, listing for each k, the 
values of co{k),... ,C7{k), and of (5.1) compared to (5.3) with D = 1000,10000, 
and 100000, and i? = 0, 1, . . . 7. 

5.1. Zeros of Pk{x). While we have managed to explicitly determine the first few 
coefficients Cr{k) of the moment polynomials Pk{x), we have not yet managed to 
understand certain aspects of these polynomials, such as the uniform asymptotics 




(5.3) 
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n 


conj. k - 


= 3 




data k = 


= 3 


conj. fc = 4 


data fc = 4 


conj. k -- 


= 5 


data A; = 5 





7.23687 X lO'-' 


7.23101 X 10" 


1.89527 X 1012 


1.88501 X 1012 


6.00428 X 


IQi^ 


5.91051 X IQi'' 


1 


1.56965 X 


10^" 


1 


57239 


X 


10^" 


5.67575 X 1012 


5.70833 X 1012 


2.45298 X 


IQi^ 


2.47886 X IQi'' 


2 


2.15687 X 


10^" 


2 


15368 


X 


IQi" 


9.17127 X 10i2 


9.12987 X 10i2 


4.68619 X 


IQi^ 


4.64908 X lOi'"^ 


3 


2.63814 X 


10^° 


2 


62463 


X 


101° 


1.24573 X lOi^ 


1.23432 X 10" 


7.10198 X 


IQi^ 


7.04187 X IQi^ 


4 


3.05562 X 


101° 


3 


06922 


X 


10i° 


1.55847 X 10" 


1.5683 X 10" 


9.63318 X 


IQi^ 


9.6445 X 10" 


5 


3.42903 X 


IQio 


Q 
O 


A A-\ AO 


X 




1.8585 X IQi^ 


1.87265 X 10" 


1.22457 X 


IQi^ 


1.24349 X 1016 


6 


3.76958 X 


lOio 


3 


76835 


X 


lOi" 


2.14798 X 1013 


2.15861 X lOi-'' 


1.4919 X 


lOi^ 


1.51619 X 1016 


/ 


4.08439 X 


lOio 


4 


05663 


X 


lOio 


o /IOC/1 \^ 1 nl3 
Z.4zo4o X iU 


o Q^om \y 1 nl3 
Z.o I zUi X iU 


L.l\)6\)o X 


1 ni6 
iU 


i.Doy 1 Z X W 


8 


4.37832 X 


IQio 


4 


39075 


X 


lOio 


2.70108 X 1013 


2.724 X 10" 


2.03988 X 


1016 


2.06017 X 1016 


9 


4.65486 X 


lOi" 


4 


65312 


X 


lOi" 


2.96679 X 10" 


2.94271 X 10" 


2.3189 X 


101^5 


2.26023 X 1016 


10 


4.91663 X 


lOio 


4 


91363 


X 


lOlo 


3.22631 X 10" 


3.24807 X 10" 


2.60051 X 


1016 


2.69184 X 1016 


11 


5.16565 X 


lOio 


5 


17448 


X 


lOio 


3.48022 X 10" 


3.47606 X 10" 


2.88433 X 


1016 


2.87018 X 1016 


12 


5.40352 X 


lOi" 


5 


39624 


X 


lOio 


3.72905 X 10" 


3.73482 X 10" 


3.17002 X 


IQifi 


3.18035 X 1016 


13 


5.63152 X 


lOlo 


5 


65418 


X 


lOio 


3.97319 X 10" 


4.00187 X 10" 


3.45733 X 


1016 


3.48184 X 1016 


14 


5.85072 X 


lOio 


5 


83654 


X 


lOio 


4.21303 X 10" 


4.1917 X 10" 


3.74603 X 


1016 


3.70813 X 1016 


15 


6.062 X 10^" 


6 


08708 


X 


lOi" 


4.44887 X 10" 


4.48257 X 10" 


4.03594 X 


1016 


4.08236 X 1016 


16 


6.2661 X 10^° 


6 


27652 


X 


lOlo 


4.68097 X 10" 


4.69566 X 10" 


4.32693 X 


1016 


4.3287 X 10" 



Table 1. This table compares the conjectured value (5.2) to ac- 
tual data (5.1) for intervals [50000n, 50000(n + l)], n = 0, 1, . . . , 16, 
and A; = 3, 4, 5. The fit is to two or three decimal places, consistent 
with the remainder stated in (1.1). 



n 


conj. k = 6 


data k = 6 


conj. k = 7 


data k = 


= 7 





2.15456 X 10" 


2.08527 X 10" 


8.45652 X 10" 


7.99015 X 


10" 


1 


1.18835 X 10" 


1.20686 X 10" 


6.24627 X 102" 


6.3773 X 


102" 


2 


2.69034 X 10" 


2.66481 X 10" 


1.67709 X 1021 


1.66563 X 


1021 


3 


4.56155 X 10" 


4.56713 X 10" 


3.18661 X 1021 


3.25679 X 


1021 


4 


6.72399 X 10" 


6.61933 X 10" 


5.1125 X 1021 


4.87831 X 


1021 


5 


9.12928 X 10" 


9.3828 X 10" 


7.42365 X 1021 


7.74635 X 


1021 


6 


1.17439 X 10" 


1.21474 X 10" 


1.00952 X 1022 


1.06992 X 


1022 


7 


1.45431 X 10" 


1.31266 X 10" 


1.31065 X 1022 


1.11053 X 


1022 


8 


1.75076 X 10" 


1.75386 X 10" 


1.64403 X 1022 


1.61306 X 


1022 


9 


2.06221 X 10" 


1.95439 X 10" 


2.00815 X 1022 


1.83038 X 


1022 


10 


2.3874 X 10" 


2.61353 X 10" 


2.40171 X 1022 


2.8627 X 


1022 


11 


2.72527 X 10" 


2.70986 X 10" 


2.82354 X 1022 


2.82074 X 


1022 


12 


3.07492 X 10" 


3.06639 X 10" 


3.2726 X 1022 


3.20372 X 


1022 


13 


3.43557 X 10" 


3.43848 X 10" 


3.74797 X 1022 


3.70176 X 


1022 


14 


3.80656 X 10" 


3.7414 X 10" 


4.24878 X 1022 


4.13975 X 


1022 


15 


4.18729 X 10" 


4.25286 X 10" 


4.77427 X 1022 


4.86676 X 


1022 


16 


4.57724 X 10" 


4.53193 X 10" 


5.32373 X 1022 


5.1628 X 


1022 


17 


4.97592 X 10" 


4.98651 X 10" 


5.89648 X 1022 


6.0058 X 


1022 



Table 2. Conjecture v.s. data for fc = 6,7, same intervals as the 
previous table. 
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Figure 1. The horizontal axis is n, and the vertical axis depicts 
(conjecture — data)/data for the values in Tables 1 and 2. For 
any finite interval, as A; ^ oo, the main contribution to the 2fcth 
moment comes from the largest value of |C(l/2 + on that 

interval. This explains the feature that, for a fixed interval, the 
actual moment tends to progressively deviate from the conjectured 
value as k increases. 



R 


cfl(.5) 


(5.3), D = 1000 


(5.3), D = 10000 


(5.3), D = 100000 





1.1299287453321533 


1463.83 


17768.4 


193494. 


1 


.19628236755422853 


1523.55 


18258.1 


197413. 


2 


.03248602185728907 


1525.93 


18271.4 


197491. 


3 


-.5289095729314908 


1516.37 


18234.3 


197335. 


4 


3.2346669444094671 


1531.24 


18275.5 


197459. 


5 


-21.381296730027876 


1505.43 


18222.2 


197343. 


6 


166.38844209028643 


1559.87 


18310.8 


197488. 


7 


-1529.2695739774642 


1419.97 


18120.1 


197237. 


(5.1) 




1521.27 


18257.1 


197425. 



Table 3. The coefficients CR{k), and conjecture v.s. data for 
fc = .5 for three intervals. The bottom row gives (5.1) for the 
interval [100, L>], with D = 1000, 10000, and 100000. For each D, 
we compare this to the value of (5.3), i? = 0, 1, . . . ,7. 



of the coefficients, or uniform asymptotics of Pk{x) with x a function of k. The 
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p 

li, 


Cfl,(3.2) 


/•c q\ n 1 nnn 

^^o.oj, jy — iuuu 


t^O.oj, U — iUUUU 


/"c o\ r) innnnn 

1^0. Oj, U — iUUUUU 





.37531596173465401 x 10 


-6 


1968.83 


1.16353 X 10*^ 


2.19960 X 10« 


1 


.34462154217944847 x 10" 


-4 


40049.5 


1.65169 X 10^ 


2.41753 X 10^ 


2 


.12662390083082525 x 10" 


-2 


336190. 


9.78885 X 10^ 


1.12289 X 10^° 


3 


.23963666452208821 x 10 


-1 


1.52891 X 10^ 


3.2097 X 10^ 


2.94868 X 10^" 


4 


.2526426167678357 




4.22213 X 10*^ 


6.6336 X 10^ 


5.06417 X 10^° 


5 


1.5214668466274718 




7.72205 X 10^ 


9.6529 X 10^ 


6.47041 X 10^° 


6 


5.3060442651520751 




1.03793 X 10^ 


1.12055 X 10^ 


7.01449 X 10^" 


7 


11.121264784324178 




1.16045 X 10^ 


1.16894 X 10^ 


7.14177 X 10^° 


(5.1) 




1.15305 X 10^ 


1.16746 X lO'^ 


7.16886 X 10^" 



Table 4. Conjecture v.s. data for k = 3.2. 



R 


Cfl(1.8) 


(5.3), D = 1000 


(5.3), D = 10000 


(5.3), D = 100000 





.13885991555298723 


15298.6 


604203. 


1.58922 X 10^ 


1 


1.2590684761107478 


46198.5 


1.42746 X 10^ 


3.20931 X 10^ 


2 


2.4174835075472416 


59612.2 


1.66821 X 10^ 


3.56224 X 10^ 


3 


2.546894763686222 


62863.9 


1.70753 X 10^ 


3.60492 X 10^ 


4 


-2.21426710514627 


62199.9 


1.7021 X 10*^ 


3.60059 X 10^ 


5 


3.223904454789757 


62432.5 


1.70339 X 10^ 


3.60134 X 10^ 


6 


46.42674651960987 


63260.1 


1.70659 X 10^ 


3.60268 X 10^ 


7 


-840.1304443557953 


59448.3 


1.69608 X 10^ 


3.59953 X 10^ 


(5.1) 




61744.5 


1.70134 X 10" 


3.60129 X 10^ 



Table 5. Conjecture v.s. data for k = 1.8. For this value of k, 
and the range we examined, J? = 4 or 5 give the best approxima- 
tion. 



R 


CR{.b + l) 


(5.3), D = 1000 


(5.3), D ^ 10000 


(5.3), D = 100000 



1 
2 
3 

4 
5 
6 
7 


1.3117481341987813 + 1.211708767666727* 
-3.0693034820213132 + 2.309977688777579i 
23.861826126198446 - 5.4045694962616631i 
-111.54278536885322 - 35.79807241977336* 

828.16689710582718 + 437.514818042632z 

-5808.11341189128 - 8339.592888954564i 
15613.29091863494 + 101218.4464636376* 

188541.27977634034 - 1175857.723687032* 


-308.872 + 439.126* 
-508.454 +246.589Z 
-335.035 + 646.347* 
-285.625+ 118.290* 
-546.679 +1199.98i 
1514.537- 1342.42i 
-6736.01 + 2796.92i 
23708.43 - 2789.22Z 


-3698.00 + 2357.78* 
-4331.26 + 957.574i 

-4243.1 + 2618. 59i 
-3667.307+ 1304.32i 
-4747.02 + 3257.57i 

-738.290 - 0.1097i 
-12358.2 + 3830.64i 

24043.14+702.364* 


-34129.1 + 8908.71* 
-35042.9 - 44.9533* 
-36981.1 + 6688.15* 
-34054.1 + 3538.21* 
-37543.5 + 6786.35* 
-30111.2 + 3451. 52i 
-45377.2 + 5624.57* 
-4932.63 + 6133.48* 


(5.1) 




-340.843 + 383.859i 


-3946.25+ 1883.17* 


-35140. + 4830.47i 



Table 6. Conjecture v.s. data for /c = .5 + *. The data here is 
not as convincing as for the other values of fc, but, nonetheless, the 
early terms do give a reasonable approximation, and we believe the 
fit would improve with more substantial data. 



latter is needed, for example, to properly understand how large |C(l/2 + it)\ can 
get [FGH]. 

5.2. Plots of Cr{k). We present in Figures 2 and 3 some graphs of the coefficients 
Cr{k) with -1/2 < fc < 11/2, for r = 0, 1 . . . , 7. 




Figure 2. Graphs of Cr{k) with -1/2 < fc < 11/2, for r = . . . , 7. 
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012345 012345 

Figure 3. The first figure depicts tlie grapli of log(|co(fc)|), 
while the next seven depict log(|cr(fc)/co(A:)|), for r = I... ,7. 
The asymptotic behaviour of log(co(A;)) as A; — *■ oo is imphed 
by [CGo] [KcS] and is, to leading order, — fc^log(fc). The cusps 
occur at zeros of Cr{k), some of which are accounted for by the 
fact that, for non-negative k gZ, Pk{x) is a polynomial of degree 
fc^ so that Cr-(fc) = if r > fc^. 
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6. Remarks about other families: orthogonal and symplectic 

In this paper we have explained two approaches for obtaining the coefficients 
c,.(fc) of Pk{x). The first involves explicitly determining the residue on the r.h.s. 
of (1.2). Theorems 1.2- 1.4, and the procedure given in Section 2 describe this 
in detail. The second approach involves using the combinatorial sum (1.11), using 
small shifts, and high precision. 

The same methods can be taken for other families of L-functions, for instance in 
determining the lower order terms in the moments of L(l/2, Xd), quadratic Dirichlet 
L-functions, or of L£;(l/2, Xd), the L-functions associated to the quadratic twists 
of a given elliptic curve, to name just two examples, in both cases evaluated at the 
critical point. The former is an example of a unitary symplectic family, while the 
latter is an example of an orthogonal family [KaS] . See [CFKRS] where we discuss 
these examples in detail. As with the Rieman zeta function, conjectures are given 
for the full asymptotics of their moments, expressed in terms of multi-dimensional 
residues and also as combinatorial sums. In that paper, we used Method 2 of 
Section 4.2 for the analogous combinatorial sums to numerically compute lower 
terms for the moments, and verify the full asymptotics. 

For the elliptic curve family, the next to leading term in the asymptotics of the 
moments has been worked out explicitly, and a test has been divised to verify the 
first two terms in the asymptotics of that particular family with an application to 
estimating the number of elliptic curves of rank greater than zero [CPRW]. See 
also [BMSW] for a survey of results related to the latter question. 
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